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We propose a model-independent scaling method to study the physical properties of high- 
temperature superconductors in the normal state. We have analyzed the experimental data of 
the c-axis resistivity, the in-plane resistivity, the Hall coefficient, the magnetic susceptibility, the 
spin-lattice relaxation rate, and the thermoelectric power using this method. It is shown that all 
these physical quantities exhibit good scaling behaviors, controlled purely by the pseudogap energy 
scale in the normal state. The doping dependence of the pseudogap obtained from this scaling anal- 
ysis agrees with the experimental results of angle-resolved photoemission and other measurements. 
It sheds light on the understanding of the basic electronic structure of high- Tc oxides. 



I. INTRODUCTION 

The mechanism of high-Tc superconductivity remains 
one of the fundamental issues unsolved in condensed 
matter physics. In particular, a unified theory toward 
the understanding of the rich phase diagram of high-Tc 
superconductors (HTSC) has not been established. At 
half-filling, the parent compounds of HTSC are antifer- 
romagnetic insulators. Upon doping, the antiferromag- 
netic long-range correlation is suppressed and the high-Tc 
superconductivity develops above a critical doping level. 
At low doping but in the metallic state, a pseudogap 
phase with missing entropy or spectra is discoveredi^'^'^ 
Around the optimal doping, a strange metal phase with 
a linear in-plane resistivity and probably a quantum crit- 
ical point emerges.^ In the heavily overdoped regime, the 
conventional Landau Fermi liquid behaviors are gradu- 
ally recovered. Throughout the whole doping range, the 
electron states seem to be intrinsically inhomogeneous,® 

A must step in the understanding of the phase diagram 
of HTSC is to identify the energy scales and the control 
parameters or interactions of low-energy excitations of 
HTSC. In the normal state at low doping, the pseudo- 
gap is believed to be one of the characteristic param- 
eters of low-lying excitations. The corresponding tem- 
perature below which the pseudogap effect is observed 
is commonly used as a boundary to separate the pseu- 
dogap and the strange metal phases, although there is 
no real phase transition between these phases. However, 
due to the uncertainty in the definition of the pseudo- 
gap, different experiments have adopted different criteria 
to determine the pseudogap energy scale. This leads to 
the claim of the existence of two pseudogaps, namely, the 
upper and lower pseudogapsji In the strange metal phase, 
the quantum critical fluctuation is strong and the tem- 
perature itself may serve as a dynamic control parameter 
as suggested by the marginal Fermi liquid theory)^ 

The scaling analysis of experimental data in HTSC is 
a simple but powerful tool in elucidating the underlying 
physics without invoking a specific model. In the normal 
state, if the pseudogap is a predominant energy scale con- 
trolling low-energy excitations, then the low-temperature 



behavior of any measurement physical quantity should 
satisfy a doping-independent scaling law, although the 
analytic expression of the scaling function is unknown. 
Based on this idea, we analyzed recently the tempera- 
ture dependence of the c-axis resistivity pc of HTSC and 
found that it obeys a universal scaling law given hy^ 

Pe(r) = ^exp(^^), (1) 

where a is a doping dependent coefhcient, T is tempera- 
ture, and A is the pseudogap. As shown in Ref. Eq. 
dl]) results from the interplay between the anisotropic c- 
axis hopping integrali^iii and the dJ.■2._y■2.-\l]^e symmetry 
of the pseudogap. It agrees excellently with the exper- 
imental data of multilayer HTSC and resolved a long- 
standing puzzle regarding the physical origin of the semi- 
conductorlike temperature dependence of pc in the pseu- 
dogap phase. Furthermore, it suggests that there is only 
one energy scale controlling the low-energy excitations 
around the antinodal points and the interlayer hopping 
within each unit cell is coherent. 

However, if the scaling function is unknown, the scal- 
ing analysis becomes difficult. This has, in fact, ham- 
pered the application of the scaling analysis. A com- 
monly adopted approach in the scaling analysis in the 
literature is to assume that by normalizing both the mea- 
surement quantity F{T) and the temperature by the cor- 
responding values at a sample-dependent characteristic 
temperature T*, then all the experimental data should 
fall onto a single curve described by the scaling function 

F(T*) ^ \x* J ^ ' 

However, in real materials, this scaling analysis often fails 
since F(T) generally contains the terms which are not 
scaling invariant. For example, the impurity contribu- 
tion to the in-plane resistivity or other physical quanti- 
ties is not scaling variant. Moreover, since T* is unknown 
prior to the analysis, this formula is difficult to be im- 
plemented practically even if it is correct. Empirically, 
T* is often determined from some special features ap- 
pearing in the measurement data. For example, for the 



2 



in-plane resistivity, T* is determined from the tempera- 
ture below which the resistivity begins to deviate from 
its high-temperature linear-T behavior. However, to de- 
termine unambiguously the linear-T region is not always 
possible since in some low doping samples, the measured 
temperature may not be high enough to reach the linear- 
T regime. Furthermore, the deviation from linear to non- 
linear T is a crossover, not a phase transition, and a small 
measurement error may result in a large error in T*. 

In this paper, we propose a scaling method and apply 
it to analyze the normal state properties of HTSC. This 
method extends our previous scaling analysis of the c- 
axis resistivity to other physical quantities. It breaks the 
barrier in the use of the simple formula ([2]) and allows a 
model-independent scaling analysis to be done reliably. 
We have reanalyzed the experimental data of the c-axis 
resistivity of HTSC with this method. By comparison 
with our previous results,^ we find that this method, in- 
deed, works very well. It provides a simple but powerful 
approach for analyzing experimental results. Further- 
more, this method is model independent. It can be ap- 
plied not only to the high- Tc cuprates, but also to any 
other materials where the single parameter scaling be- 
havior is valid. 

This paper is arranged as follows. Section |TT] gives 
an introduction to the scaling method based on a least 
square fit of an unknown scaling function to the experi- 
mental data. In Sec. IIIIl we apply the method to analyze 
the scaling behaviors of a number of physical quantities of 
HTSC in the normal state, using the experimental data 
published in the literature. In Sec. IIVI we analyze the 
universal behavior of the energy scales obtained in Sec. 
mil and discuss its physical implications. Section IVl gives 
a brief summary. 



II. METHOD OF SCALING ANALYSIS 

In this section, we present a generic method for ana- 
lyzing the scaling behavior of a set of experimental data 
{Fi{T), i — 1, N}. Here, T can be temperature, pres- 
sure, external field, or any other controllable variable 
used in experiments. The subscript i is a sample index to 
which the physical quantity T as a function of T is mea- 
sured experimentally. N is the total number of samples. 
In the discussion below, in order to be directly relevant to 
the scaling analysis presented in Sec. IIIIl we assume that 
T is the temperature and a sample represents a specified 
high- Tc compound. 

We start by assuming that in a relevant temperature 
range the low-lying physics is governed only by one en- 
ergy scale A. Thus the measured physical quantity F[T) 
satisfies a simple scaling law 



(3) 



where a, (3, and A are all doping dependent, but temper- 
ature independent. For the physical quantities to be dis- 



cussed in Sec. IIIIl (3 is generally the contribution of impu- 
rities or other extrinsic interactions. A is a characteristic 
energy scale of the system. It controls the dynamics of 
the system. 3-{T) is a universal (doping-independent) 
scaling function. Its temperature dependence is deter- 
mined by the low- lying excitations and interactions. 

The scaling function !F{T) is generally unknown. This 
is the difficulty commonly met in the data analysis. How- 
ever, the scaling method introduced here does not depend 
on the detailed formula of T{T) , provided that the single- 
parameter scaling hypothesis Eq. ^ is valid. This is a 
merit of the method. It provides a simple but powerful 
approach to probe the intrinsic connection between dif- 
ferent samples and to determine the doping dependence 
of the characteristic energy scale A without invoking a 
specific model. 

The aim of the scaling analysis is to determine from 
the measurement data the scaling parameters (a^, A^, (3i) 
and the optimized scaling functions T(x) so that all the 
data can be rescaled onto a universal curve. This can, in 
principle, be achieved by minimizing the total deviation 
of the scaling function between any two samples: 



N Nk 



i<j k 



where 



[F{A,T)-Pi 



(4) 



(5) 



Fi{A.iT) is the value of F of the zth sample at tem- 
perature AiT. Nk is the number of sampled tempera- 
ture points used in optimizing the scaling function. Nk 
can be adjusted in the minimization. Initially, can 
take roughly the value of average measured temperature 
points. 

However, due to the scaling behavior of T{T), not all 
the parameters {ai, Ai, (3i) can be uniquely determined 
by the minimization of 5F if T is unknown. Indeed, 
from Eqs. @ and ([5]), it can be shown that if {cki. A,, 
(i = 1, N) is a set of parameters minimizing then 
{ai/ Us, Ai/ As, I3i — PsCHi/as) with arbitrary but nonsin- 
gular (tts. As, Ps) will also minimize ST. (as, As, /3s) are 
unknown and can be taken as the scaling parameters of a 
reference sample. This means that only the relative val- 
ues of ai, Ai, and f3i with respect to a reference sample 



A. - ^, (6) 

as 

B. = (7) 

a = p^-Aps, (8) 

can be fixed. From the definitions (O-®, it is straight- 
forward to show that As = Bs = 1 and Cs = 0. 

Using the relative scaling parameters {Ai, Bi, Ci), one 
can define a new scaling function 



/ 



asT 



T 

'ba. 



(9) 
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F{T) can then be expressed as 



F(T)-A/(|)+C. 



(10) 



For the reference sample, the scahng function is the mea- 
surement curve itself: F{T) — f{T). 

Equation PH)) is nothing but to scale all experimental 
data onto the measurement curve of the reference sample. 
This suggests that the reference sample should be chosen 
such that its temperature interval is broad enough to 
cover the whole temperature range physically interesting 
and the data quality is among one of the best. 

The relative scaling parameters {Ai^ Bi,Ci} can now 
be determined by minimizing the total deviation of the 
scaling functions: 



where 



i<j k 



MT) = - [F,(B,T) 



(11) 



(12) 



In general, BiT may not be exactly the temperature point 
experimentally measured. The value of Fi(BiT) can be 
obtained from the measurement data by interpolation, 
provided that BiT is within the measured temperature 
interval. For each pair of fi and fj, if BiT^ or BjTk 
is outside the measured temperature interval for the zth 
or jth sample, then the corresponding term in Eq. Q 
should be excluded from the summation. The minimiza- 
tion of Sf can be done, for example, using the standard 
subroutine given in Ref. IT^ . 



TABLE I: The scaling parameters for the in-plane resistivity 
pab of YBa2Cu306+4 published by Ito et al}^ The parameters 
of the reference sample, S = 0.68, are in bold face. 



Sample (5) 


A 
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C (IQ-'^flcm) 


0.68 
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FIG. 1: (a) The experimental data of the in-plane resistivit 
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pabiT) of YBa2Cu306+i , extracted from Fig. 2(a) of Ref . 
(b) The scaling plot of the experimental data. The reference 
sample is (5 = 0.68. 



Below we take the in-plane resistivity data of 
YBa2Cu306+5 published by Ito et al}^ as an example 
to show how the method works. For simplicity, here we 
only use the experimental data for three of the samples, 
5 = (0.68,0.78,0.85). Figure [2a) shows the measure- 
ment data of Pab{T) for these three samples. In the scal- 
ing analysis, we take 5 = 0.68 as the reference sample 
and Nk to be roughly equal to the measured tempera- 
ture points of the reference sample. The relative scaling 
parameters can then be determined by numerically min- 
imizing Eq. (|lip . The resulting scaling curves are shown 
in Fig. [ijb) and the scaling parameters are given in Table 
III 



In Eq. (fTU|) . if C ^ F{T), then the scaling equation is 
approximately given by 

F{T) « Af . (13) 

In this case, the scaling equation can be reexpressed as 



F{T) ^ f{T/B) _ 
F{T') f{T*/B) 




where T' is a sample-dependent characteristic tempera- 
ture. g{T) is a rescaled function of f{T). Equation 
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is precisely the scaling equation defined by Eq. It is 
a commonly used scaling equation in the scaling analy- 
sis of experimental data. However, it should be pointed 
out that this equation is valid only when the temperature 
independent term C can be safely ignored or reliably sub- 
tracted from F{T) in Eq. 



III. SCALING ANALYSIS OF EXPERIMENTAL 
DATA 

In this section, we apply the scaling method to an- 
alyze the experimental data of high- Tc cuprates in the 
normal state. We will first analyze the scaling behav- 
ior of the c-axis resistivity pc- Since an approximate but 
accurate expression for the scaling function of pc is avail- 
able, this allows us to determine the absolute values of 
the scaling parameters. For other measurement quan- 
tities, including the in-plane resistivity pab{T), the Hall 
coefficient Rh{T), the magnetic susceptibility x{T), the 
spin-lattice relaxation rate l/TiT, and the thermoelec- 
tric power S{T), only the relative scaling parameters can 
be determined. 

Our scaling analysis is based on the experimental data 
already published in the literature. We collect as much as 
we can the experimental data of HTSC from which a sys- 
tematical analysis of the scaling behaviors can be done. 
The preference is given to the latest published data if 
there are considerable differences between the data pub- 
lished by different groups. The chemical formula of the 
compounds with their abbreviations analyzed in this pa- 
per are given in Table |TT1 The scaling analysis here will 
be limited to the superconducting samples in the normal 
state. The data for the nonsuperconducting samples will 
not be analyzed. In all figures and tables presented in 
this paper, the parameters for the reference samples will 
be in bold face to distinguish them from other parame- 
ters. 

TABLE II: HTSC compounds and their abbreviations ana- 



lyzed in this paper. 

YBaaCusOs+a Y123 

Yo.8Cao.2Ba2Cu306+5 Ca-Y123 

YBa2Cu408 Y124 

BiaSraCaCuaOg+i Bi2212 

BiaSraCaaCuaOio+a Bi2223 

La2-2:Sr2-a:Cu04 La214 

BiaSra-^La^CuOs+i La-Bi2201 

HgBa2Cu04 Hgl201 

TlSraCaCuaOr-i T11212 



In the comparison of the scaling parameters for differ- 
ent families of HTSC, we will use the superconducting 
transition temperature Tc and its empirical formula pro- 



posed by Presland et al.^ 

-^-l-82.6(p-0.16f, (15) 

c.max 

to determine the effective carrier concentration p. Here, 
Tc,max is the maximal superconducting transition tem- 
perature. For La2-2;Sr2;Cu04, the carrier concentration 
is equal to the doping concentration of Sr ions, p ~ x. 

A. c-axis resistivity Pc{T) 

In the pseudogap phase of high- Tc cuprates, the c-axis 
resistivity pc behaves very differently from its in-plane 
counterpart pab- Along the Cu02 plane, pab shows a 
metal-like temperature dependence. It decreases with 
decreasing temperature. However, along the c-axis, pc 
behaves as a semiconductor. It increases with decreasing 
temperature. 

This dramatic difference between pc and pab is not 
what one might expect within conventional Fermi liq- 
uid theory. To resolve this issue, a number of theoreti- 
cal models based on the dynamic confinement of charge- 
spin separated particle a^^'^^i"'^'^ or the incoherent inter- 
layer hoppingi^ii^ were proposed. Most of the theories 
predicted that pc should diverge in a certain power law 
of T at low temperature. However, it seems that none of 
these theories can account quantitatively or even quali- 
tatively the experimental data. 

The semiconductorlike behavior of pc results, as we 
recently pointed out;^ from the interplay between the 
da.2_j^2-like pseudogap and the anisotropic c-axis hopping 
integraljiSiii p^ contributes mainly from the quasiparti- 
cles around the antinodal points. The nodal contribution 
is completely suppressed by the interlayer hopping matrix 
elements. Since the pseudogap is a prevailing energy scale 
governing the c-axis dynamics in the pseudogap phase, it 
is natural to assume that the c-axis resistivity satisfies a 
scaling law governed purely by the pseudogap A. 

In cuprate superconductors, if Cu atoms in the two 
neighboring Cu02 planes lie collinearly along the c-axis, 
then the interlayer hopping integral between these two 
planes is given by 

tc ^ {coskx — cos ky )^, (16) 

where (kx, ky) are the in-plane momenta of electrons, tc 
vanishes along the nodal direction. Based on this for- 
mula, we showed in Ref. ^ that for multilayer cuprates, 
Pc is approximately given by 

Pc{T) = aTc (^1) , (17) 

where 

Tc{x) = xexp [ — ) . (18) 
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FIG. 2: The scaling function Tc of the c-axis resistivity Pc{T) 
for Y123, Bi2212, and Bi2223, compared with the theoretical 
curve determined by Eq. (|18|l (dashed line). The experimen- 
tal data are obtained from Refs. 22 23 24 for Bi2212, from 
Ref.liifor Bi2223, and from Refs. l2fl and ^211 for Y123. Bi2223 
samples are labelled with the notations defined in Ref. [25l . 

Equation. (|17p is a special case of Eq. It holds when 
the residual resistivity contributed by disorder scattering 
f3 is vanishingly small compared with the contribution of 
pseudogap to pc- 

The scaling function (fT8|) . as shown in Fig. [21 agrees 
excellently with the measurement data for Y123 pub- 
lished by Yan et al3 and by Babic et al.^^ for Bi2212 by 
Watanabe et aA^^^^^ and for Bi2223 by Fujii et al^ 

The values of the scaling parameters a and A are given 
in Ref. The pseudogap A, as shown in Fig. [3l 

drops almost linearly with doping. This doping depen- 
dence of the pseudogap agrees with the angle-resolved 
photoemission^ (ARPES) as well as other measurement 
data.— >^ The values of A (not shown here) for two over- 
doped samples of Bi2212 from Chen et al— deviate obvi- 
ously from the other points. However, their pc data can 
be well scaled onto the universal curve, as shown in Fig. 
3(a) of Ref. [oj. A probable explanation is that the true 
doping levels for these two samples may not be as high 
as reported. The four samples, F, G, H and I, of Bi2223 
near optimal doping (here, we adopt the notations used 
in Ref. I25i) have almost the same Tc- However, their A 
obtained from pc are very different. This might be due 
to the inhomogeneity of charge carriers in these samples, 
since Tc is determined mainly by the fraction of a sample 
where the carrier concentration is close to the optimal 
doping, while pc is the contribution of the whole sample. 

The existence of the universal scaling law of pc, es- 



FIG. 3: The doping dependence of the pseudogap A obtained 
from the scaling analysis of the experimental data shown in 

tj2]with Eqs. ([T3 and mil. [a], [b], and [c] refer to Refs. 
I23I and [24 . and i2Q, respectively. 

pecially its activated behavior, implies that the c-axis 
hopping is predominantly coherent, rather than incoher- 
ent as usually believed. The reason is actually simple. 
If the interlayer hopping is incoherent, then the excita- 
tions around the gap nodes may have substantial con- 
tribution to Pc which may break this scaling law. For 
multiple layer cuprates, the intralayer coupling may be 
different from the interlayer coupling. However, the dif- 
ferent coupling between Cu02 planes does not change the 
fact that the pseudogap is the only energy scale govern- 
ing the quasiparticle excitations around the antinodes in 
the normal state. Therefore, Eq. (fT7|) holds irrespective 
of the number of Cu02 planes in each unit cell. 

TABLE III: The scaling parameters, Ac and Be, of the c- 
axis resistivity pc defined by Eq. (f20)l for La-Bi2201. The 
experimental data and the carrier concentration p are taken 
from Ref. 30. 



Sample (x) 


Tc (K) 


Doping (p) 


Ac 


Be 


0.84 


1.4 


0.10 


2.50 


1.13 


0.73 


14 


0.11 


1.37 


1.14 


0.66 


23 


0.12 


1.0 


1.0 


0.49 


31 


0.14 


0.79 


0.51 


0.39 


38 


0.16 


0.67 


0.41 



The above scaling analysis is done based on Eq. (flSl) 
by assuming that the interlayer hopping integral is given 
by Eq. (|16p . However, in single-layer cuprate com- 
pounds, Cu atoms of two adjacent Cu02 planes do not 
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FIG. 4: Scaling behavior of the c-axis resistivity pc{T) for the 
single-layer cuprate superconductor La-Bi2201. The experi- 
mental data were taken from Ref.jsO- The dashed line denotes 
the scaling function defined by Eq. ([18}. 



lie collinearly along the c-axis. In this case, the c-axis 
hopping integral becomesii 



tc ^ (cos kx — cos ky)"^ cos cos 



(19) 



2 2 

It vanishes along both the nodal and antinodal direc- 



tions. The c-axis hopping, thus, contributes mainly from 
the low-lying excitations between the nodal and antin- 
odal points. As the pseudogap is the dominant energy 
scale in the pseudogap phase, the scaling law of the c- 
axis resistivity [Eq. ([T7])] should still hold. However, the 
scaling function of the single-layer cuprates is expected 
to be different. In this case, an accurate expression for 
the scaling function is not available. Thus, we are unable 
to determine the absolute values of the scaling param- 
eters A and a from the scaling analysis. However, the 
relative scaling parameters can be determined using the 
scaling method introduced in the previous section. 

We have analyzed the scaling behavior of pc for La- 
Bi2201 (Ref. HO) using the formula 



^cfc 



T 



(20) 



Again, the contribution from impurity scattering to pc 
is ignored since it is much smaller than the contribution 
from the pseudogap effect. 

Figure. [4] shows the scaling function fc for La-Bi2201. 
The scaling parameters are given in Table [1111 The analy- 
sis shows that pc, indeed, exhibits a good scaling behavior 
in this one-layer material. However, the scaling function 
is different from that for the multilayer materials at low 
temperatures. It is also different from the logarithmic 
divergence behavior as observed by Ando et al. in La214 
cuprates 1^ 

The scaling behavior of pc for both the single- and 
multiple-layer cuprates indicates that the interlayer dy- 
namics is, indeed, governed by the interplay between the 
anisotropic interlayer hopping integral and the pseudo- 
gap effect in all high-Tc cuprates. 



J 



B. In-plane resistivity pab{T) 

In contrast to the semiconductor like behavior of pc, the 
in-plane resistivity pab of high-Tc cuprate superconduc- 
tors is metal-like in the normal state. In the underdoped 
and optimally doped materials, Pab exhibits a univer- 
sal linear behavior at high temperatures The phonon 
scattering can lead to a linear resistivity. However, the 
Debye temperature determined by applying the Bloch- 
Gruneisen formula is too low too account for the experi- 
mental data of Bi2201.'^"' The linear behavior of pab could 
be a manifestation of strong correlations. It is a char- 
acteristic behavior of marginal Fermi liquid,^' where the 
inelastic scattering rate scales linearly with temperature. 
It may also result from gauge^^i^i^i or quantum critical 
fluctuations. 

Below a characteristic temperature T*, pab begins to 



deviate from the linear behavior. This deviation is corre- 
lated with the pseudogap elfect, and T* is believed to be 
the onset temperature below which the pseudogap opens. 
However, the opening of the pseudogap does not lead to 
a phase transition. There are no thermal anomalies ob- 
served in the specific heat or other thermodynamic quan- 
tities around T* . 

To analyze the scaling behavior of the in-plane resis- 
tivity, we assume pab to satisfy the following scaling law: 

Pab{T) = Aabfab ( ^) + Cab- (21) 
\^ab / 

Here, the residual resistivity Cab should be retained since 
the impurity contribution to pab is no longer negligible 
compared with the inelastic contribution of electrons to 

Pab- 

We have applied Eq. ([^ to the experimental data 
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FIG. 5: The scaling function of the in-plane resistivity for ¥123" (Ref.[Tl), Y123'' (Ref.lS^ 
and La214 (Ref. 1381). The Bi2223 samples are labeled using the notations given in Ref. |25 



Bi2212 (Ref.lla), Bi2223 (Ref. 25), 



published by Ito et alJ^ for Y123, by Watanabe et ali2^ 
for Bi2212, by Fujii et al.^^ for Bi2223, and by Nakano 
et al^ for La214. Figure [5] shows the scaling function 
for these compounds. The corresponding scaling param- 
eters are shown in Table ITVl The measurement data (not 
shown in the figure) begin to deviate from the universal 
scaling curves near Tc due to superconducting fluctua- 
tions. 

The result of Fig. [5] shows that pab, indeed, satisfies 
the simple scaling law described by Eq. (HH). Moreover, 
all the curves shown in Fig. [5] are obtained by taking 
Y123 6 = 0.68 as a reference. This means that pab can 
be scaled onto a single curve for all these materials. Thus 
the scaling function of pab is universal. It does not de- 
pend on the chemical structure nor on the doping level. 



This suggests that the in-plane resistivity is governed by 
the same scattering mechanism in all cuprate supercon- 
ductors 

The striking scaling behavior of pab indicates that the 
characteristic temperature T* above which pab varies lin- 
early with temperature is proportional to the scaling pa- 
rameter Bab- In Ref. [H, we showed that T* is propor- 
tional to the pseudogap A determined from the c-axis re- 
sistivity within measurement errors, independent of dop- 
ing. Thus, Bab is also proportional to A. This means 
that, same as for the c-axis resistivity, the pseudogap A 
is also a control energy scale for the in-plane resistivity, 
although Pab is mainly the contribution of nodal quasi- 
particle excitations. 

Wuyts et ali^i^ did a similar scaling analysis for the 
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TABLE IV: The scaling parameters, Aab, Bab, and Cab for 
the in-plane resistivity of Y123°_(Ref. Y123'' (Ref;_p2), 
Bi2212 (Ref. Hi), Bi2223 (Ref. [2|), and La214 (Ref. [s^as 
shown in Fig. \5\ The notations used in Ref. '25' are adopted 
to label the samples of Bi2223. Y123° with 5 = 0.68 is taken 
as the reference sample. The unit of Cab is 10~'*f2cm. The 
doping p of ¥123", Bi2212 and Bi2223 are obtained from the 
empirical formula Eq. (|15|) . with the maximum Tc being 93.54 
K for (5 = 0.90 of ¥123", 89.0 K for Bi2212 at 6 = 0.22, and 
108.0 K for F sample of Bi2223. The values of for La214, 
¥123" Bi2212, and Bi2223 are obtained from Ref. 0, [H, HI, 
andi25l. respectively. 



In the overdoped regime, pab is not linear-T dependent 
in nearly the whole temperature range. "^^ The experimen- 
tal data do not fall onto the scale curves as shown in 
Fig. [5l This change of the temperature behavior of pat in 
the overdoped regime can be understood from the change 
of the Fermi surface topology revealed by ARPES.*- In 
the overdoped region, the Fermi surface becomes elec- 
tronlike. This may affect strongly the dynamic behavior 
of electrons in the Cu02 planes. 



¥123" 



¥123" 



Bi2212 



Bi2223 



La214 



Sample 


Tc (K) 


Doping 


(p) Aab 


Bab 


Cab 


0.45 


54.96 


0.088 


3.83 


1.91 


0.09 


0.58 


64.67 


0.098 


1.48 


1.22 


0.13 


0.68 


67.04 


0.100 


1.0 


1.0 


0.0 


0.78 


80.26 


0.118 


0.63 


0.81 


0.0074 


0.85 


92.08 


0.146 


0.39 


0.60 


-0.10 


0.90 


93.54 


0.16 


0.22 


0.50 


-0.19 


0.55 






2.02 


1.33 


0.19 


0.60 






1.57 


1.20 


-0.31 


0.65 






1.35 


1.15 


-0.22 


0.70 






0.89 


0.92 


-0.095 


0.75 






0.70 


0.81 


-0.070 


0.80 






0.55 


0.65 


-0.26 


0.85 






0.52 


0.78 


0.049 


0.2135 


71.0 


0.11 


1.18 


0.96 


1.70 


0.217 


77.0 


0.119 


1.07 


0.89 


1.10 


0.22 


83.0 


0.131 


0.65 


0.73 


1.12 


0.24 


89.0 


0.16 


0.51 


0.63 


0.60 


0.245 


87.86 


0.173 


0.48 


0.82 


0.76 


0.255 


87.4 


0.175 


0.32 


0.54 


0.24 


B 


93.0 


0.118 


0.12 


1.09 


0.31 


D 


104.0 


0.139 


0.10 


0.99 


-0.20 


E 


106.0 


0.145 


0.087 


0.81 


-0.69 


F 


108.0 


0.16 


0.07 


0.69 


-0.90 


0.10 


30.75 


0.10 


0.84 


2.06 


9.47 


0.125 


32.0 


0.125 


0.60 


1.97 


6.93 


0.14 


36.62 


0.14 


0.51 


2.05 


5.77 



in-plane resistivity of Y123. However, the characteris- 
tic temperatures (or energy scales) they determined are 
not very accurate. Their scaling curves of pab (Fig- 1 of 
Ref. and Fig. 7 of Ref. |40|) do not look as good as 
those shown in Fig. [5l 

In Fig. [HKd) , we only show the experimental data for 
slightly underdoped La214 samples (x = 0.1,0.125, and 
0.14). For heavily underdoped samples {x = 0.04,0.06, 
and 0.08), we find that the data deviate significantly from 
the universal scaling curve below T*. This deviation was 
observed only in the La214 samples. It may be due to 
the suppression of the scattering rate by the formation of 
stripe or other competing orders in these compounds4i 
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FIG. 6: The scaling functions of the Hall coefficient Rh{T) 
for (a)La214 (Ref SI), (b) La-Bi2201 (ReLlH) and (c) ¥123 
(Ref. 1481 ). The dashed line in (a) is obtained from the scaling 
curves determined by Hwang et al. (Ref. .50. ) . 
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C. Hall coefficient Rh{T) 

The Hall coefficient Rh{T) is an important quantity 
in characterizing the nature of charge carriers. In a con- 
ventional metal with spheric Fermi surface and isotropic 
scattering rate, the Hall coefficient is inversely propor- 
tional to the carrier concentration, independent of tem- 
perature. The sign of Rh reflects the type of conducting 
charge carriers. Rh is negative or positive if the charge 
carriers are electrons or holes. However, in doped tran- 
sition metal oxides, such as high- Tc cuprates, the Hall 
coefficient is strongly temperature dependent.— "^i^ It is 
determined not just by the carrier concentration, but also 
by the scattering rates and the curvature of the Fermi 
surface. Other effects, such as magnetic skew scattering, 
can also affect the temperature dependence of Rh^ 

In HTSC, Rh shows a strong temperature and doping 
dependence. At high temperature, Rh increases rapidly 
with decreasing temperature. After reaching a maxi- 
mum, Rh drops down to low temperature in most of the 
samples. This is the typical temperature dependence of 
Rh in HTSC. It was observed in Y123,i^ ha2U^^^ 
Bi systems^i^i^ Hgl212,55 and Tl systems-^iSi^ 

The complex temperature dependence of the Hall co- 
efficient remains one of the hardest problems to be re- 
solved. Within the theory of charge-spin separation, 
AndersonSi^ proposed to use the Hall angle Qh in- 
stead of the Hall coefficient Rh to understand the Hall 
anomaly. He argued that due to the charge-spin sep- 
aration, the Hall angle, which is defined by the ratio 
between the transverse and longitudinal conductivities, 
Qh = axy/cTxx, should be determined purely by 

the transverse relaxation rate (i.e., the relaxation rate 
perpendicular to the Fermi surface). This eliminates the 
ambiguity in the explanation of the Hall coefficient, since 
it is determined by both the longitudinal and transverse 
relaxation rates. Anderson further argued that as the 
transverse relaxation rate is determined by the spin exci- 
tations, which is relatively normal, the Hall angle should 
follow the temperature dependence of normal Fermi liq- 
uid, i.e.. 



cote 



H 



aT^ + 7, 



(22) 



where a is a temperature independent coefficient and 7 is 
the impurity contribution.^ » This quadratic temperature 
dependence of the Hall angle, indeed, agrees with the 
experimental observation at optimal doping. However, 
in both underdoped and overdoped regimeS )^^'^°i^^'^^i^^ 
the temperature exponent deviates generally from 2, and 
the above expression breaks down. 

To reveal the physics behind the anomalous temper- 
ature dependence of the Hall effect without invoking a 
specific model, we have analyzed the scaling behavior of 
Rh- We assume the scaling function of Rh to have the 
form defined as in Eq. ^T^ 



TABLE V: The scaling parameters for the Hall coefficient 
Rh shown in Fig. [g] The unit of Ch is 10"^ cm^C"^ Tc for 
La214 are obtained by interpolating the data given in Ref. 0. 
The hole concentration p for La-Bi2201 and Y123 are obtained 
from Eq. psp with the maximum Tc,max ~ 33.03 K for x — 
0.44 of La-Bi2201, and 93 K for 5 = 0.95 of Y123. 





Sample 


Tc (K) 


Doping 


: (p) Ah 


Bh 


Ch 


La214 


0.21 


30.99 


0.21 


0.10 


0.52 


0.34 




0.18 


37.14 


0.18 


0.25 


1.08 


0.23 




0.15 


37.93 


0.15 


0.54 


1.40 


0.035 




0.12 


31.88 


0.12 


0.96 


1.62 


-0.42 




0.08 


22.3 


0.08 


2.23 


2.15 


-2.04 




0.05 


2.66 


0.05 


4.34 


2.68 


-5.62 


La-Bi2201 


0.24 


24.06 


0.218 


0.14 


0.91 


0.91 






ou.u 




n 99 


n 87 

U.O ( 


111 

J Li. 




0.44 


33.03 


0.16 


0.35 


1.12 


1.00 




0.57 


28.44 


0.118 


0.49 


1.23 


1.08 




0.66 


19.95 


0.09 


0.61 


1.27 


1.70 


Y123 


0.55 


55.0 


0.088 


1.76 


1.27 


0.56 




0.6 


57.0 


0.091 


1.52 


1.15 


0.20 




0.65 


58 


0.092 


1.46 


1.09 


-0.035 




o.r 


60.0 


0.094 


1.0 


1.0 


0.0 




0.75 


62.0 


0.096 


0.86 


1.05 


-0.20 




0.8 


69.0 


0.103 


0.67 


1.08 


-0.11 




0.85 


83.0 


0.123 


0.70 


0.96 


-0.17 




0.95 


93.0 


0.16 


0.35 


0.78 


0.073 




7.0 


91.0 


0.176 


0.28 


0.67 


0.031 



Rh{T) = AhJh 



T 



C 



H- 



(23) 



By fitting the experimental data with this formula using 
the method introduced in Sec. [TTl the scahng parameters 
(Ah, Bh,Ch) and the scaling function fn can then be 
determined. 

Figure [6] shows the scaling function fn for Y123, 
La214, and La-Bi2201. The corresponding scaling pa- 
rameters are given in Table [V] The experimentaj_data 
were extracted from Ref. for Y123, from Ref. |4^ for 
La214, and from Ref . [sl for La-Bi2201. The experimen- 
tal data of Y123 with S < 0.55 are not included since 
the temperature range measured is too narrow to allow 
a reliable scaling analysis to be done. 

For La214, we find that Rh exhibits a good scaling be- 
havior in nearly the whole temperature range, as shown 
in Fig.[ni;a). This is also true for La-Bi2201 [Fig.li^b)]. It 
suggests that the dynamical behavior of Rh is still gov- 
erned by a single energy scale. However, in contrast to 
the in-plane resistivity, the scaling curve of La214 cannot 
be perfectly scaled onto the scaling curve of La-Bi2201, 
except in an intermediate temperature regime. For Y123, 
the experimental data can also be scaled onto a single 
curve in the high-temperature regime, above the peak 
temperature of Rh- However, at low temperatures, fn 
shows very different temperature dependence for different 
dopings. This difference might be caused by the contri- 
bution of CuO chains in Y123. Electrons in CuO chains 
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are more disordered in the underdoped samples than in 
the optimally doped one. 

A similar scaling analysis has been done by Hwang 
et alM- and by Chen et al^ The scaling equation they 
used is essentially the same as Eq. (1^51) , but with different 
notations. In their notation, the scaling equation is given 
by 

RH{T)^R^ + R*Hf(^^y (24) 

is the high-temperature value of Rh, which is ap- 
proximately temperature independent at high tempera- 
tures. T* is a characteristic temperature to be deter- 
mined. Empirically, they assumed T* to be the crossover 
temperature from a temperature dependent to a tem- 
perature independent Rh at high temperature. How- 
ever, the crossover temperature (if exists) is very high, 
well above the temperature range they measured, in the 
underdoped samples. They cannot determine reliably 
the crossover temperature, even by extrapolation. Thus, 



D. Thermoelectric power S{T) 

The thermoelectric power or the Seebeck coefhcient 
S(T) is one of the transport quantities complementary 
to the resistivity and Hall effect. It reveals the proper- 
ties of quasiparticle excitations both near and away from 
the Fermi level. It can be used to judge whether the 
charge carriers are electrons or holes from the sign of 
S{T). It can also be used to quantify the charge carrier 
concentration. In HTSC, empirically, the value of S{T) 
at T = 290 K was found to be a good measure of the hole 
concentration!^ 

In high- Tc oxides, the thermopower S{T) is small in 
the superconducting state due to the suppression of the 
pairing gap to the quasiparticle excitations In the nor- 
mal state, S{T) increases with temperature. It exhibits 
a maximum and then drops monotonically at high tem- 
peratures. At a given temperature, S{T) decreases with 
increasing doping. For most of the high-Tc compounds, 
including Bi22 12, §2*64^^66 La-Bi2201,^ TIUU^^ 
and Hgl201,^?^i^ S{T) is positive in the underdoped 
regime, but becomes negative in the overdoped regime. 
At optimal doping, S{T) becomes negative at high tem- 
perature. However, in La214, S{T) was found to be pos- 
itive in the whole doping range . ^^i^^'^'^ The thermopower 
of overdoped Y-based cuprates also behaves differently 
from other compounds. It shows a positive slope, in con- 
trast to the negative slope in other compounds, at high 
temperature. This is probably due to the contribution of 
CuO chains. 

We have applied our scaling approach to the thermo- 



their scaling analysis cannot be applied to the under- 
doped samples. This is not a problem in our approach. 
In Fig. El^a), we compared the scaling curve obtained by 
Hwang et al^ with ours. We find that these two scal- 
ing curves agree well with each other for La214 above 
the peak temperature of Rh- This suggests that the 
crossover temperature they determined is proportional 
to the parameter Bh as we determined here. 

Recently, Gor'kov and Teitel'baum.^^ and Ono et al^ 
analyzed the high-temperature behavior of Rh using a 
two-band model in La214. They assumed the high- 
temperature data of Rh to be thermally activated, re- 
sulting from strong charge fluctuations between the ef- 
fective lower and upper Hubbard bands. They found 
that the high-temperature data of Rh can indeed by ex- 
plained by this simple picture. Their results suggest that 
the thermal excitation gap between the lower and upper 
Hubbard bands is significantly smaller than the (direct) 
optical charge transfer gap. 



I 

electric power. The scaling formula is given by 



S{T) = Asfs [^]+Cs. (25) 



Figure. [7] shows the scaling curves for Y123, Ca-Y123, 
Hgl201, and La214. The scaling parameters are listed 
in Tables IVII and IVIII The experimental data were ex- 
tracted from Ref. i72| for Y123, from Ref. (zlfor Ca-Y123, 
from Ref. [6^ for Hgl201, and from Refs. l73l and [t^ for 
La214. The scaling curves for two heavily underdoped 
samples with Tc < 2 K and three heavily overdoped sam- 
ples of Hgl201 (Ref. [6^ deviate significantly from the 
universal scaling curve and are not included in the figure. 
For Y123 and Ca-Y123, we only analyze the underdoped 
samples since in the overdoped regime the chain contri- 
bution becomes important, which breaks the scaling law. 
The chain contribution can, in fact, be seen already in the 
shghtly underdoped sample of Y123 [5 ~ 0.175), whose 
high-temperature data of S{T) already begin to deviate 
away from the scaling curve at high temperatures. For 
some samples of Hgl201 and Ca-Y123, the measurement 
data fall faster than the universal scaling curves at low 
temperatures. This can be attributed to the supercon- 
ducting fluctuations. 

The scaling behavior of S{T) in HTSC has been ex- 
tensively studied by a number of groups using the scal- 
ing formula like that defined in Eq. ^,66 ,69^,72 Qui- 
scaling curves are consistent with their results. However, 
our data are much less scattered than theirs. 
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FIG. 7: The scaling function of the thermoelectric power S{T) for Y123 (Ref.lzi), Ca-Y123 (Ref.lzi), Hgl201 (Ref.lH), La214" 
(Ref.lli), and La214'' (Ref.[72). The notations defined in Ref . |6i3 are used to label the Hgl201 samples. Sample 8 of Hgl201 
is the reference. 



E. Uniform magnetic susceptibility x{T) 

The uniform magnetic susceptibility x(T) measures 
basically the density of states at the Fermi level in con- 
ventional Landau Fermi liquid. However, in high- Tc cop- 
per oxides, the susceptibility is strongly affected by an- 
tiferromgantic spin fluctuations. The parent compounds 
of HTSC are half- filled antiferromagnetic Mott insulators 
with long-range Neel order Upon doping, the Neel or- 
der is rapidly suppressed, but antiferromagnetic fluctua- 
tions persist up to slightly overdoping.i^LZ^ 

At half- filling, x{T) shows a sharp peak around the 
Neel temperaturei ^^i^° This peak shifts down to lower 



temperature with doping and disappears completely 
when the superconductivity emerges. In the normal 
state, the magnetic susceptibility x(^) first increases 
with increasing temperature, develops a broad peak, and 
then drops down at high temperaturei^i2^i^>22i^i2ii2^ 

In real materials, the magnetic susceptibility is 
strongly affected by magnetic impurities. These impu- 
rities contribute a Curie term to Xj which diverges as 
1/T at low temperatures. The impurity contribution is 
strongly sample dependent. In order to analyze the in- 
trinsic behavior of the magnetic susceptibility, this Curie 
term of impurities should be subtracted from the raw 
data first. 
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TABLE VI: The scaling parameters for the thermoelectric 
power S{T) of Y123, Ca-Y123, and Hgl201 in Fig. [3 The 



unit of Cs 


is fiV/K 


. Sample 


8 of Hf 


;1201 


is taken as 


a refer- 


encG. 


















Sample 


T. (K) 


Doping 


(P) 


As 


Bs 


Cs 


Y123 


0.65 








75.53 


1.80 


4.17 




0.60 


13.31 


0.057 




60.37 


1.60 


3.22 




0.53 


44.66 


0.079 




42.42 


1.27 


2.46 




0.44 


55.51 


0.089 




15.67 


0.99 


2.36 




0.35 


59.8 


0.093 




10.36 


0.94 


2.03 




0.24 


74.4 


0.109 




1.18 


0.83 


1.46 




0.175 


90.7 


0.14 




-1.46 


0.79 


0.48 


Ca-Y123 


0.38 


85.50 


0.16 




0.31 


0.61 


-1.17 




0.43 


85.00 


0.151 




0.41 


0.80 


-1.07 




0.51 


81.58 


0.136 




0.99 


0.81 


-1.57 




0.54 


78.55 


0.128 




1.10 


0.84 


1.13 




0.57 


77.19 


0.125 




1.15 


0.85 


4.21 




0.69 


47.04 


0.085 




1.46 


1.05 


12.86 




0.99 


37.69 


0.077 




2.01 


1.40 


30.07 


Hgl201: 


3 


26.0 


0.05 




3.65 


1.58 


33.56 




4 


46.0 


0.057 




3.81 


1.43 


11.79 




5 


62.0 


0.069 




2.50 


1.33 


15.24 




6 


72.0 


0.09 




1.98 


1.26 


4.88 




7 


77.0 


0.103 




1.90 


1.12 


-3.27 




8 


83.0 


0.11 




1.0 


1.0 


0.0 




9 


91.0 


0.119 




1.18 


1.07 


-4.67 




10 


95.0 


0.127 




0.96 


1.01 


-2.80 




11 


98.0 


0.157 




0.68 


0.98 


-3.01 



To elucidate the intrinsic property of the magnetic sus- 
ceptibihty, we have analyzed the scaling behavior of x(T) 
with the following single-parameter scaling equation: 



T 



(26) 



Figure [8] shows the scaling curves of x(^) for La214 
and Bi2212. The experimental data were extracted from 
those published by Nakano et al.^^ for La214 and by 
Watanabe et alM- for Bi2212. The corresponding scal- 
ing parameters are shown in Table IVIIII In the scal- 
ing analysis for the La214 samples with x = 0.20, 0.22, 
and 0.26, a Curie term C/T is subtracted from the ex- 
perimental data; the corresponding values of C, i.e., 
C = 12.41,37.23, and 60.23 (in units of lO"'^ emu/g), 
were obtained by Nakano et al^ 

We find that the susceptibility for both La214 and 
Bi2212 exhibits a good scaling behavior. For Bi2212, 
x(T) in heavily overdoped samples begins to deviate from 
the scaling curve near Tc. This is likely to be due to 
strong superconducting fluctuations. 

Our universal scaling curve is consistent with the scal- 
ing analysis given by Johnston^^ and Nakano et al.'^^ for 
La214, and by AUgeier and Schillings^ for Bi2212. The 



TABLE VII: The scaling parameters for the thermoelectric 
power ^(r) of La214 in Fig. [T] The data are referenced to 
sample 8 of Hgl201 shown in Table ED The unit of Cs is 
/iV/K. Tc of La214'' is obtained from the data published in 
Ref. by interpolation. Tc of La214'' is obtained from Fig. 1 
of Ref. III. 





Sample 


Tc (K) 


Doping 


:(p) As 


Bs 


Cs 


La214'' 


U, 


1 c 

.io 


37.93 


0.15 


1.23 


0.87 


16.44 




U, 


.io 


37.14 


0.18 


0.90 


0.69 


8.07 




0, 


.2 


34.62 


0.2 


0.74 


0.63 


5.70 




U, 


.21 


30.99 


0.21 


0.64 


0.60 


4.46 




U, 


.26 


23.77 


0.23 


0.52 


0.61 


1.95 




0, 


.25 


17.29 


0.25 


0.45 


0.65 


0.15 




0, 


.26 


14.4 


0.26 


0.42 


0.72 


-1.48 




0, 


.28 


8.08 


0.28 


0.44 


0.67 


-0.85 


La214'' 


0, 


.06 


5.0 


0.06 


53.51 


2.01 


4.18 




0, 


.08 


22.3 


0.08 


51.78 


1.74 


3.22 




0, 


.1 


30.75 


0.1 


40.03 


1.54 


2.53 




0, 


.125 


32.0 


0.125 


31.55 


1.17 


1.66 




0, 


.135 


35.9 


0.135 


30.04 


1.01 


0.88 




0, 


.15 


37.93 


0.15 


18.66 


0.78 


0.79 




0, 


.17 


38.4 


0.17 


17.73 


0.69 


0.87 




0, 


.2 


34.62 


0.2 


10.59 


0.62 


0.83 




0, 


.22 


27.81 


0.22 


6.34 


0.54 


0.61 




0, 


.24 


19.96 


0.24 


3.41 


0.54 


0.53 




0, 


.27 


12.0 


0.27 


2.69 


0.53 


0.43 



scaling analysis of Nakano et al. was made based on 
the scaling formula defined by Eq. ([14]). ''^ The problem 
with that kind of analysis is that the characteristic tem- 
perature T* defined in Eq. has to be determined 
empirically prior to the scaling analysis. The character- 
istic temperature T* for the susceptibility was generally 
determined from the peak temperature of x- However, in 
heavily overdoped samples, no peak structure has been 
observed within the whole temperature measured. This 
has limited the application of that kind of scaling analy- 
sis. In addition, to fully satisfy Eq. a constant term 
also needs to be subtracted for each set of data. This 
is also difficult if the measured temperature range is not 
broad enough. Nevertheless, we find that their scaling 
curves, as shown in Fig. [SI agree well with ours. 

The scaling analyses given by Johnston^ and by AU- 
geier and Schilling^^ are based on the high-temperature 
series expansion for a two-dimensional antiferromag- 
netic Heisenberg model. The scaling function obtained 
by Johnston deviates slightly from the universal scal- 
ing curve obtained at high temperature. The high- 
temperature scaling curve, as shown by Johnston^^ and 
by AUgeier and Schilling, agrees with the temperature 
dependence of the susceptibility of the two-dimensional 
antiferromagnetic Heisenberg model without doping. 

In obtaining the scaling function, a temperature inde- 
pendent term is subtracted from xiT). This term can be 
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FIG. 9: The temperature independent magnetic susceptibility 
Xo ~ Cx + ^xXo.s where xo.s ~ 1-4 x 10~^emu/g for La214 
with X — 0.18 as obtained from Nakano et al. (Ref. 38i). 



dau diamagnetism, and the Pauli paramagnetism of the 
band electrons if the effect of antiferromagnetic correla- 
tions is ignored. The core diamagnetism and the Van 
Vleck paramagnetism are doping independent. The Lan- 
dau diamagnetism is generally small. Thus, the doping 
dependence of xo is mainly affected by the Pauli suscep- 
tibility, which is proportional to the density of states at 
the Fermi level. Hence, the change of xo with doping 
will correspond to the change of the density of states at 
the Fermi level. This simple observation is consistent 
with the result of Allgeier and Schilling^ as well as the 
measurement of the specific heati^ 



FIG. 8: The scaling function of the uniform magnetic suscep- 
tibility x(r) for La214 (Ref. (H) and Bi2212 (Ref. HI). The 
solid, dash-dot and dash lines represent the rescaled sca ling 
curves obtained by Johnston fRef. |83) ■ Nakano et al fRef. l38f) 
and Allgeier and Schilling (Ref. ISa). respectively. 

expressed as 

Xo=C^+A^Xo,s, (27) 

where Xo,s is the value of xo for the reference sample. The 
variation of xo with doping may reflect the change of the 
density of states at the Fermi level. Thus, it is interesting 
to analyze the doping dependence of this term. 

Figure [9] shows the doping dependence of xo ■ Xo in- 
creases monotonically with doping. This temperature in- 
dependent term might be the contribution of the core 
diamagnetism, the Van Vleck paramagnetism, the Lan- 



F. Spin-lattice relatxation rate 1/Ti 

The nuclear magnetic resonance probes the local spin 
dynamics via the measurement of the Knight shift, the 
spin-lattice relaxation rate 1 /Ti , and other spin response 
functions.*® The Knight shift measures the shift of the 
resonance frequency induced by the conduction electrons. 
It is proportional to the uniform magnetic susceptibility 
x(T). The temperature dependence of the Knight shift 
should follow the scaling law of the uniform susceptibility 
as discussed in the previous section. This is, indeed, sup- 
ported by the experimental measurement (see, for exam- 
ple, Ref. 'sT). Below we will discuss the scaling behavior 
of the spin-lattice relaxation rate 1 /Ti . 

The high-temperature dependence of l/TiT of ®^Cu 
shows a Curie- Weiss-like behavior i^^i*^ This can be at- 
tributed to the contribution of antiferromagnetic fluctu- 
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TABLE VIII: The scaling parameters for the magnetic sus- 
ceptibihty x{T) shown in Fig. [S] The unit of is 10"^ 
emu/g. Tc of La214 is obtained from Ref. by interpola- 
tion. In obtaining the scaling curves, a Curie term C/T is 
subtracted from the measurement data for the La214 samples 
with X = 0.20, 0.22, and 0.26. The corresponding values of C, 
obtained from Fig. 8 of Ref. H, are 12.41, 37.23, and 60.23 
lO"'^ emu/g, respectively. 





Sample 


Tc (K) 


Doping (p) 






Cx 


La214 


0.08 


22.3 


0.08 


1.40 


5.97 


-1.60 




0.10 


30.75 


0.10 


1.41 


4.18 


-1.57 




0.14 


36.62 


0.14 


1.28 


2.31 


-1.11 




0.18 


37.45 


0.18 


1.0 


1.0 


0.0 




0.20 


34.62 


0.20 


1.04 


0.68 


0.31 




0.22 


27.81 


0.22 


1.24 


0.53 


-0.001 




0.26 


14.57 


0.26 


1.38 


0.46 


-0.027 


Bi2212 


0.22 


82.0 


0.13 


1.33 


3.04 


-1.29 




0.25 


89.57 


0.16 


1.10 


2.22 


-0.66 




0.26 


87.19 


0.18 


0.66 


1.60 


0.61 




0.27 


83.45 


0.186 


0.72 


1.18 


0.62 




0.28 


78.7 


0.195 


0.51 


0.74 


1.42 




0.30 


68.51 


0.22 


0.66 


0.40 


1.84 



ations since 1/TiT of ^^Cu is dominated by spin fluctu- 
ations near Q = (tTjTt). It exhibits a broad maximum in 
an intermediate temperature regime and then drops at 
low temperature. 

TABLE IX: The fitting parameters. An, Bti, and Cti for 
the scaling analysis of 1/T\T shown in Fig. 1101 The unit of 
Cti for ''^'Cu of Bi2212 and T11212 is s"^K-\ and for ^''Y of 

1 

S 





Sample 


Tc (K) Doping 


(p) Ati 


Bti 


Cti 


Bi2212 ^^'Cu 


0.125 


79.0 


0.13 


0.86 


1.31 


0.72 




0.20 


86.0 


0.16 


1.0 


1.0 


0.0 




0.225 


77.3 


0.20 


1.17 


0.85 


-1.67 


T11212 ^^Cu 




70.0 




2.44 


0.72 


-4.72 






54.0 




2.29 


0.56 


-3.05 






10.0 




0.48 


0.16 


16.18 


Y123 


1.0 


88.8 


0.19 


0.075 


0.88 


1.51 




l-£ 






0.069 


1.35 


1.53 




0.85 


90.3 


0.14 


0.085 


1.63 


1.15 




0.75 


67.9 


0.10 


0.087 


1.50 


0.83 




0.63 


57.3 


0.09 


0.088 


2.30 


0.90 




0.53 


52.5 


0.086 


0.097 


2.80 


0.80 




0.48 


38.0 


0.07 


0.070 


3.35 


0.77 




0.41 


15.0 


0.06 


0.028 


3.74 


0.68 



We have analyzed the scaling behavior of the exper- 
imental data of 1/TiT published by Ishida et al.^'^ for 
63Cu in Bi2212, by Magishi et alS^ for ^^Cn in T11212, 
and by AUoul et ali^ for in Y123. The scaling equa- 
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FIG. 10: The scaling function for l/TiT of '*='Cu in Bi2212 
and T11212, and that of *^Y in Y123. The experimental data 
are obtained from Ref. 190 for Bi2212, from Ref. iH for T11212, 
and from Ref. [l| for Y123. The notations of the samples are 
the same as those in the corresponding references. The Bi2212 
OP86 is taken as the reference sample. 



tion is assumed to be 



1 , , f T 
— — — AtiJti — 

111 \DT1 



Cti- 



(28) 



Figure [TO] shows the temperature dependence of the scal- 
ing function fTi- The corresponding scaling parameters 
are given in Table. IIXI 

Within experimental errors, we find that 1/TiT shows 
a good scaling behavior. For all the materials shown in 
Fig. [TOl the data of 1/TiT can be scaled on a common 
curve in a relatively wide range of temperature. This, 
again, suggests that the normal state dynamics is con- 
trolled by a single energy scale. This single-parameter 
scaling behavior is consistent with the existing theory 
of spin fluctuations, such as the antiferromagnetic Fermi 
liquid theory proposed by Millis et al^ 



IV. ANALYSIS OF THE CHARACTERISTIC 
ENERGY SCALE 

In the preceding section, we have analyzed the scaling 
behaviors of the transport coefficients, including the re- 
sistivity, the Hall effect, and the thermoelectric power, 
and the magnetic response functions, including the spin 
susceptibility and the spin-lattice relaxation rate, in the 
normal state of HTSC. These coefficients probe differ- 
ent aspects of low-energy excitations and are physically 
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distinct. However, we find that they aU show good scal- 
ing behaviors. For most of the measurement quantities, 
the temperature dependence of the corresponding scahng 
functions are universal, depending neither on the doping 
concentration nor on the chemical structure of the mate- 
rials measured. 

Among the three scaling parameters, the energy scale 
A or the relative energy scale B defined in Eq. ([7]) is 
the most important one. It characterizes the basic en- 
ergy scale governing the temperature dependence of a 
response function. For the c-axis resistivity, we have de- 
termined the absolute values of A using the approximate 
scaling function of pc derived in our previous work for 
multilayer cuprates.~ For other measurement quantities, 
as the analytic formula of the scaling functions are un- 
known, only the ratio of A with respect to a reference 
sample, B — A/A^, is determined. Nevertheless, we find 
that these characteristic energy parameters determined 
from different coefficients show a common trend with 
doping. As shown in Tables IIIHIVIII they all decrease 
with increasing doping. 



900 



600 



300 




FIG. 11: The doping dependence of the energy scales obtained 
with Eq. (|29[) . The values of the scaling factor riy are listed 
in Table [X] The dashed line is the linear fit to the scaled 
pseudogap energies obtained by various experimental probes. 
Here, we take the pseudogap energy obtained from the c-axis 
resistivity of Bi2212 as a reference, as shown in Fig. |3l 

The similar doping dependence of the characteristic en- 
ergy parameters suggests that these relative energy scales 
obtained by different probes may have a common phys- 
ical origin. This can be examined by rescaling all the 
relative energy scales with respect to the absolute energy 
scale of the pseudogap A obtained from the scaling anal- 
ysis of Pc- They should fall onto a single curve if they 
are, indeed, the pseudogap energies. To do this, let us 



TABLE X: The values of the scaling factor r/y for different 
probes y = pc, Pab,RH,X^ l/TiT, and S of the cuprates. The 
superscripts a and 6 for S and Rh of La214 refer to the data 
published in Ref. fiol and Ref. fsS, respectively. 





Pc 


Pab 


Rh 


S 


X 




Bi2212 


1 


455 






125 


267 


Bi2223 


1.04 


375 










Y123 


1.67 


465 


450 


450 




182 


Ca-Y123 








435 






La214 




405 


225" 


255" 


150 










0.46'' 


315'' 






La-Bi2201 


450 




310 









Hgl201 



385 



introduce the following formula: 



(29) 



where rjy is a scaling factor and Ay is the characteristic 
energy scale. The subscript y represents the measured 
physical quantity, i.e., y = Pc, Pab, Rh, S, x, and l/TiT. 
The scaling factors rjy can be determined by the least 
square fit using the approach introduced in Sec. |TT1 

The fitting parameters of rjy are given in Table 1x1 By 
substituting them into Eq. (pQ]) . we can obtain the val- 
ues of Ay. The result, as shown in Fig. [TTl indicates 
that all the energy scales determined from the scaling 
analysis given in Sec. IIIII have the same doping depen- 
dence within experimental errors, which result mainly 
from the uncertainty in the determination of doping con- 
centration. It suggests that all the dynamic coefficients 
analyzed in Sec. lIIII are. indeed, governed by the same en- 
ergy scale. This is a remarkable result since different co- 
efficients probe different responses of charge and/or spin 
degrees of freedom. For example, the c-axis resistivity p^ 
is susceptive to the charged excitations around the antin- 
odal points, while the in-plane resistivity pab is mainly af- 
fected by the scattering of charged quasiparticles around 
the node points. The uniform magnetic susceptibility 
probes the spin fluctuations around k = (0,0), while the 
spin-lattice relaxation rate is strongly dependent on the 
antiferromagnetic fluctuations around k = (7r,7r). 

In the underdoped regime, the characteristic energy 
A drops almost linearly with doping. This is consistent 
with the doping dependence of the pseudogap observed 
by ARPES)^ tunneling)2^ and other measurementsi^ 
Thus, the control energy scale in this regime is, indeed, 
the pseudogap. If we extrapolate the underdoped data 
of A to zero doping, we find that A is in order of J, and 
to higher doping, A vanishes roughly at p ^ 0.23-0.25. 

In the literature, two kinds of pseudogaps (or the on- 
set temperatures of pseudogap) , which were often quoted 
as the "large" and "small" pseudogaps, were reportedii 
The large pseudogap generally refers to the characteris- 
tic temperature T* measured, for example, by the mag- 
netic susceptibility's, the Hall coefficienti^ The small 
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one could be the energy scale probed by other exper- 
imental techniques, such as the transport properties 
measurements^ or the leading-edge shift around the 
antinodal direction measured by the ARPES.— Our scal- 
ing analysis indicates that these two energy scales are, in 
fact, physically indistinguishable. 

In the overdoped regime, the universal scaling behavior 
generally breaks down. For overdoped La214, Ay deter- 
mined from the thermopower and the Hall coefficients 
separates into two branches with different energies. Ay 
determined from the thermopower is higher than that 
from the Hall effect. This may imply the existence of 
two energy scales in the overdoped regime. However, as 
the single-parameter scaling law still holds in this regime 
for both the thermopower and the Hall coefficients, fur- 
ther investigation on this issue is desired. 

Recently, there is a surge of interest in the discussion 
of two-gap energies, namely, the pseudogap in the nor- 
mal state and the relatively smaller superconducting gap 
in the superconducting state. The existence of a distinct 
energy scale in the superconducting state, whose doping 
dependence is different from that of the pseudogap, was 
first reported in the penetration depth measurements^^ 
and later in the Andreev reflection measurements This 
lower energy scale characterizes the low-lying excitations 
around the gap nodes and appears only in the supercon- 
ducting state. Recently, ARPES,^^ Raman scattering,— 
and inelastic neutron measurements^^ have further con- 
firmed the existence of these two distinguished energy 
scales. Furthermore, the ARPES has revealed that the 
two-gap structure is intimately connected with the arc 
(or pocket) feature of the Fermi surface of HTSC in the 
underdoped regime. It is believed that the superconduct- 
ing gap develops predominately on the Fermi arc below 
Tc- This is consistent with the early ARPES measure- 
ment datai2S However, there exists also other experimen- 
tal measurements, which suggest that there is only one 
energy scale in the superconducting state and the super- 
conducting gap around the nodal points is nothing but 
an extension of the pseudogap in the arc area in the su- 
perconducting statCf^ This one energy scale scenario 
is consistent with the picture of resonant valence bonds 



based on the charge-spin separationiSi as well as that of 
preformed pairs. ^''^ 

In this work, we have shown that there is only one 
energy scale in the normal state. However, as we have 
only analyzed the scaling behavior of the experimental 
data in the normal state, we are unable to address the 
issue of two energy gaps in the superconducting state. It 
is of great interest to extend the single-parameter scal- 
ing method introduced in Sec.|TT]to the superconducting 
state at which two energy scales (or control parameters) 
may exist. This would then allow us to judge whether 
there is only one or two energy scales in the supercon- 
ducting state from the model-independent scaling analy- 
sis of various transport and thermodynamic coefficients. 
V. SUMMARY 

We have introduced a scaling method to study the scal- 
ing behavior in the normal state of HTSC. We have an- 
alyzed the scaling behavior of the c-axis resistivity, the 
in-plane resistivity, the Hall coefficient, the thermoelec- 
tric power, the magnetic susceptibility, and the nuclear 
magnetic resonance, and extracted the corresponding en- 
ergy scales. It is found that all these quantities, no mat- 
ter how different they are, exhibit universal scaling be- 
haviors, controlled by a single energy scale in the nor- 
mal state. Furthermore, we find that all these energy 
scales obtained from different physical coefficients have 
the same doping dependence as the pseudogap. It shows 
that the pseudogap is the only characteristic energy gov- 
erning the low-lying excitations in the normal state of 
HTSC. 

The scaling method we introduced in Sec. [IT] is model 
independent. It provides a simple but powerful tool to 
analyze the scaling behavior of experimental data. It can 
be applied not only to the normal state of HTSC, but also 
to any other systems where the single-parameter scaling 
hypothesis, i.e. Eq. ^ or (fTU)l. is valid. 
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